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Abstract 

Using an effective lagrangian description we analyze possible new physics contri- 
butions to the most relevant muon number violating processes: fi — > e 7 and fi-e 
conversion in nuclei. We identify a general class of models in which those processes 
are generated at one loop level and in which fi-e conversion is enhanced with re- 
spect to fi — > e 7 by a large ln(m^/A 2 ), where A is the scale responsible for the new 
physics. For this wide class of models bounds on fx— e conversion constrain the scale 
of new physics more stringently than fi — ► e 7 already at present and, with the ex- 
pected improvements in fi-e conversion experiments, will push it upwards by about 
one order of magnitude more. To illustrate this general result we give an explicit 
model containing a doubly charged scalar and derive new bounds on its couplings 
to the leptons. 



1 Introduction 



The precision reached in the last years in the experiments searching for fi—e 
conversion in nuclei at PSI [1,2] and TRIUMF [3] and the expected improve- 
ment in the sensitivity of the experiments at PSI in the next years by more 
than two orders of magnitude [4] will make fi-e conversion the main test of 
muon flavour conservation for most of the extensions of the standard model 
(SM). Moreover, according to the recent BNL proposal [5] further improve- 
ments in the experimental sensitivity down to the level 10 -16 are feasible. 
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There are three processes for which there are very good experimental bounds: 
/j, — > e 7, /i-e conversion and fi — > eee. None of them occurs in the SM without 
extending it with right-handed neutrinos or extra scalars [6] but in general 
they appear in any extension of the SM in which lepton flavour conservation 
is not imposed by hand, /i — > e 7, having the photon on mass-shell, can only 
be originated from an off-diagonal (in generation space) magnetic moment. 
This interaction can only appear, in any renormalizable model, at the loop 
level. On the contrary, /i-e conversion and ji — > eee can be generated at three 
level from renormalizable interactions by exchange of scalars or gauge bosons. 
However, in many models those couplings conserve muon flavour or do not 
appear in both lepton and quark sectors. In this case /i-e conversion and/or 
/i — > eee are also generated at the loop level. If this is the case, and given 
the present experimental accuracy, one often finds [7] that the bounds on new 
physics coming from /i — > e 7 are stronger than the bounds found from /i-e 
conversion or /i — > eee. 

However, this does not need to be always the case. It could happen that the 
form factors contributing to /i-e conversion are enhanced with respect to the 
ones contributing to /i — > e 7. Indeed, it has been noticed already in the early 
works of ref. [8] that in some cases /i-e conversion constrains new physics 
more stringently than /i — > e 7. In this Letter we investigate the conditions 
in which this happens in the framework of the effective quantum field theory 
which allows one to classify in a simple way the contributions coming from a 
large variety of extensions of the SM. We point out a general class of models 
in which /i-e conversion is enhanced by large logarithms. As an example we 
will study a simple extension of the SM with a doubly charged scalar singlet 
coupled to the right-handed leptons and derive new limits on its couplings. The 
same bounds, although derived for a singlet, hold with a good precision also 
for models containing scalar triplets which are introduced usually to generate 
Majorana neutrino masses [9] and appear naturally in left-right models [10]. 
Finally we comment on other presently popular models (e.g., with broken 
-R-parity [11] or leptoquarks [12]) with the similar feature. 



2 Effective lagrangian description of theories 



Assuming that the relevant physics responsible for muon-number nonconser- 
vation occurs at some scale A > = Fermi scale, we can write the relevant 
effective lagrangian at low energies as 

c eff = c L + c R + C aL + C aR + C LL + C RR + ---, (1) 
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where 



a R 

C R = j^f^e e^ lv e jR d,F^ , (3) 

£ aL = J^^ e ^<V* P^lF^ + h.c. , (4) 
a aR 

^ L = -^(^ L )(^e iL ), (6) 

^ R = ^(e^e c fcfi )(^). (7) 



Here, as in the rest of the paper, we will assume that repeated indices, Lorentz, 
fj,, z/, or generation indices, i,j,k,l, are summed. When possible we will use 
also matrix notation in generation space, and em are chiral charged-lepton 
fields, e c L R = (eL t n) c are the charge conjugated fields and p =0 + ie /L. The 
four-fermion couplings cxf^j and a R k R j are symmetric with respect to the ex- 
changes i <-> k and/or I <-> j. 

We expect that the terms C L , C R , L aL , are generated at one loop in the 
renormalizable theories since they cannot be obtained from renormalizable 
vertices at tree level, that is the reason we already included a factor (4-7r) 2 
in the denominator. C L and Cf L will arise, for instance, in models with an 
extra scalar triplet with hypercharge 1 [9,13] or an extra scalar singlet [14,15] 
coupled to the leptonic doublet. C R and C aR will arise, for instance, in models 
with a doubly charged scalar singlet [16] coupled to the singlet right-handed 
leptons, we will study this model more carefully latter on. 

Note the particular form we have written the magnetic-moment type operators 
involving only chiral fields. The two operators C aL and C aR could be combined 
by using the equations of motion for the lepton fields. We have 

% pe L = M e e R , i pe R = M\e L , (8) 

where M e is the charged lepton mass matrix and we have used a matrix nota- 
tion to suppress generation indices. By using the equations of motion, which is 
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perfectly allowed in an effective lagrangian [f 7] at the lowest order, we obtain 
C° = C° L + C° R = (4 ^ 2A2 e eza^ (a aL M e + M e a° R ) e R + h.c. (9) 

In fact for applications we will use C a written in this form. Note that we could 
use from the beginning this form for C a as the starting effective lagrangian, but 
this is not completely equivalent to what we did. By doing that we would have 
no reason to choose the particular form of magnetic moments proportional to 
the fermion masses present in eq. (9). However, in chiral theories, like the ones 
we want to consider, magnetic moments appear always in the form eq. (9) and 
are proportional to the fermion masses. In more general theories with chirality 
explicitly broken independently of the fermion masses, operators like eq. (9) 
but with an arbitrary matrix could arise. 

Note also the form in which we have written the four fermion operators in 
terms of conjugate fields. Those operators are equivalent, after a Fierz trans- 
formation, to the usual vector-vector four fermion operators, for instance 

(eifle£ fl ) (e^e jR ) = ^(e^^e jR ) (e^7 M e iR ) . (10) 

We have chosen this form because it is simpler for loop calculations since it 
leads to only one penguin diagram while the vector- vector interaction leads to 
two types of penguin diagrams. Moreover, these operators arise naturally in 
the class of models we will consider. 

Four-charged-fermion interactions violating generation-number conservation 
however can be generated easily at tree level in a large class of models (ex- 
cluding supersymmetry with conserved R-parity). For instance C LL will appear 
in scalar models in which the scalars couple to the lepton doublet and models 
with a scalar triplet with hypercharge 1. Notice that a singly charged singlet 
cannot generate these kind of couplings, it only generates couplings with two 
charged leptons and two neutrinos [14,15]. C RR will appear in models with a 
doubly charged singlet (more on this later). Four-fermion couplings involving 
both, charged leptons and neutrinos have not been included because, as we 
will see later, they do not lead to any logarithmic enhancement of the /i-e 
conversion rate. Four-fermion couplings involving leptons and quarks could 
also be included, if not bounded already by another reasons, and would lead 
to a similar logarithmic enhancement to the one we are going to study. 

In addition to the couplings we have considered there could be four fermion 
operators involving both left-handed and right-handed fields. They could be 
generated, for instance, by exchange of Higgs doublets. However, these cou- 
plings are usually suppressed by the masses of the fermions. The consequences 
of these kind of operators at tree level have already been considered in [18]. 
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Therefore, for simplicity, we are not going to consider them in this paper. 
Moreover, a direct Z^el^^ej vertex can in principle be generated at tree level 
in some models in which ordinary fermions mix with other fermions with ex- 
otic hypercharges and at one loop in models with non-decoupling physics in 
the same way that extra couplings Z ' y&irf^hi, arise in the SM when one tries 
to make the top-quark mass very heavy [19]. In particular those couplings will 
arise when a Dirac mass term for the neutrinos is made large [20]. In these 
kind of models /i-e conversion could be sizable with respect to /i — > e 7 studied 
in ref. [21]. Since this mechanism has already been studied elsewhere [20,22] 
we are not going to consider it here anymore and will concentrate in models 
in which /i-e conversion proceeds through the photonic mechanism, that is, 
by exchange of a photon between the leptonic and the hadronic currents. 

It is important to note that the lagrangian (1) has to be interpreted as a la- 
grangian in the effective field theory approach^]. This means that four-fermion 
interactions, which are generated at tree level, can be used at one loop and 
will generate non-analytical contributions to the electromagnetic form factors. 
In fact, as we will see, those non-analytical contributions are quite indepen- 
dent of the model and are the key of the possible enhancement of the form 
factors contributing to /i-e conversion with respect to those contributing to 
/t — > e 7. If there are logarithmic contributions to the /i-e conversion rate, 
they will dominate, and since they can be computed in the effective theory, 
they are quite independent on the details of the full theory from which the 
effective lagrangian is originated. In section 4 we show, by using an explicit 
model, how this works and how the /i— e conversion rate is quite independent 
on the details of the model. 



3 /t-e conversion versus /i — > e7 

Theory of /i-e conversion in nuclei was first studied by Weinberg and Feinberg 
in ref. [23]. Since then various nuclear models and approximations are used 
in the literature to calculate coherent /i-e conversion nuclear form factors. It 
is important to note that the results from the shell model [24], local density 
approximation [25] as well as the quasi-particle RPA approximation [26] do 
not differ significantly from each other for both ^Ti and g2 8 -P& nuclei showing 
consistency in the understanding of the nuclear physics involved [26] . We follow 
the notation of ref. [25] and take into account corrections to the local density 
approximation from the exact calculations performed in the same work. The 
corrections are negligible for as the local density approximation works 
better for light nuclei but are sizable for ]^Pb. 

3 For an detailed explanation of this approach in the case of a model with a singly 
charged scalar see for instance [15]. 
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The relevant [i-e conversion matrix element can be expressed as 

47ra (1) Gp (2) 



where q is the momentum transfer and in a good approximation q 2 m —mft- 
The first term in eq. (11) describes photonic and the second term non-photonic 
conversion mechanisms, jh 2 n and ^ represent the leptonic and hadronic 
currents, respectively. The non-photonic mechanism is mediated by heavy par- 
ticles and, therefore, suppressed compared with the photonic mechanism. The 
non-photonic mode is of interest if the conversion can occur at tree level like 
in models with non-diagonal Z' [22] and Higgs [18] couplings, models with 
broken i?-parity [27], models with leptoquarks [28] or if the loop contributions 
are enhanced by some other mechanism, e.g., models with non-decoupling of 
massive neutrinos [20]. Since we do not consider tree level lepton number vi- 
olation via four fermion operators involving quarks in this Letter we shall 
concentrate in the following on the photonic mechanism only. 

Generally, the leptonic current for the photonic mechanism can be parametri- 
zed as 



J(l)=u(p e ) 



[JEO + 75/M0) lv\g ~ 



+ (^1+75^1)^^ — 



(12) 



where p e , p^ are the lepton momenta and the form factors /^ , /ei, /mo an d 
/mi can be computed from the underlying theory. The coherent /i-e conversion 
branching ratio R^ e can be expressed as [25] 



R 



fie 



c 



PeE ( 



\F( Pe 



capt 



?0 ' 



(13) 



where E e is the electron energy, T capt is the total muon capture rate, |F(p e )| 2 
is the nuclear matrix element squared and 

fo = \fE0 + /Mil' + |/ S 1 + /M0| 2 (14) 



shows the /x-e conversion dependence on the form factors. The expression for 
\F(p e )\ 2 in the local density approximation, the correction factors C to the 
approximation (compared with the exact calculation) as well as all numerical 
values of the above defined quantities for f 2 Ti and l^Pb can be found in 
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refs. [25,29]. The result reads 



8a 5 mlZ* f Z\FJp e )\ 2 & , x 



where C Ti = 1.0, C Pb = 1.4, = 17.61, = 33.81, r^ pt = 2.59 • 

10 6 s _1 , = 1.3 • 10 7 s _1 and the proton nuclear form factors are F p Tl (q) = 

0.55 and F p (g) = 0.25. Presently SINDRUM II experiment is running on 
gold [30], 7q 9 Au, but gold is not explicitly treated in ref. [25]. However, since 
Z Au = 79 and Z Pb = 82 are so close to each other then, within errors, all the 
needed quantities for H 9 Au and f^Pb are approximately equal 4 . This result is 
strongly supported by theoretical calculations and experimental measurements 
of the total muon capture rate of Pb and Au [29] . In the following we use the 
same experimental and theoretical input for both }l 9 Au and f^Pb. 

One should note that the \i — > e 7 branching ratio, 

^ = S$(i/Mir+i^r), (i6) 

depends on a different combination of the form factors. 




Fig. 1. The one-loop contribution to the electromagnetic current in the effective 
theory. 

To show the power of the effective lagrangian description of new physics we 
compute the form factors starting from the lagrangian (1). The one loop di- 
agram giving rise to /i-e conversion is depicted in fig. 1. By using the MS 
renormalization scheme and choosing the renormalization scale, /i — A, where 
A is the scale of new physics, we find for the form factors: 



Ieo — 



-r 



(4tt) 2 2A 2 



<LL , RR 



We thank H.C. Chiang and E. Oset for clarifying us this point. 
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MO 



/mi 
Jei 



(4tt) 2 2A 2 



+ « - o 



(4vr) 2 A 2 



(47r) 2 A 2 V eM eM / ' 



(17) 

(18) 
(19) 
(20) 



where r fc = m|/(— g 2 ), with m fc being the masses of the fermions running in 
the loop, and 



F(r) = lnr + 4r h (1 - 2r)Vl + 4r In . r 

3 \ yl + 4r — 1/ 



(21) 



There are three important limiting cases. If r — > (i.e., k = e) then F(r) = 
—5/3, if r ps 1 (i.e., k = /i) then F(r) ~ 0.18 and if r ^> 1 (i.e., k — r) then 
F(r) = lnr. The loop diagram in fig. 1 give contributions only to the form 
factors /eo and /mo but not to /ei and /mi- Because of the UV divergence 
find in the loop calculation those contributions always contain a term which is 
proportional to ln(g 2 /A 2 ) or ln(m 2 /A 2 ). This term which is completely inde- 
pendent of the details of the model that originate the four-fermion interaction 
gives a large enhancement for the form factors /^o and /mo while the en- 
hancement is absent in the form factors $ei and /mi- Consequently, the /i-e 
conversion is enhanced while \x — > e 7 is not. In this class of models, in which 
H~e conversion is dominated by this large logarithmic term one can neglect 
all the non- logarithmic contributions which are the ones that depend on the 
details of the complete theory. 



4 An explicit model 



Let us consider for a moment an extension of the SM by adding just a doubly 
charged scalar singlet k ++ . Its coupling to right-handed leptons are described 
by 

C-k = hije^ejR k ++ + h.c. (22) 
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Here the Yukawa coupling matrix hij is symmetric in the generation indices 
From this interaction we obtain easily, see fig. 2, the four-fermion interaction 

^Ki h iA^ e kR)(^R e jR) ■ ( 23 ) 




Fig. 2. Tree-level matching. 

This interaction is of the type C RR . Comparing with eq. (7) one can immedi- 
ately identify 



a?kfij = h*khij and A = m K . (24) 



On the other hand, by using the techniques in [15] one can easily obtain the 
contributions from matching to the full theory at one loop to the rest of the en's. 
The one loop diagrams involving k ++ are depicted in fig. 3. Those diagrams 
are computed by using dimensional regularization and, after subtraction of 
the effective theory contributions, fig. 1, they can be expanded in 1/m 2 .. We 
keep at most terms of order 1/m 2 .. At this order there are contributions to the 
self-energies and to the vertex of the photon. Those give three types of opera- 
tors: charge radius operators, eq. (2) and eq. (3), magnetic moment operators, 
eq. (4) and eq. (5), and operators that involve three covariant derivatives of 
the fermions. The last operators can be removed by using the equations of 
motion in favour of mass terms and do not lead to any interesting physics. 
Therefore, after wave function renormalization, in order to write the kinetic 
terms in canonical form, the only operators generated in this model are those 
appearing in the lagrangian (1), but only the right-handed components. By 
using the MS renormalization scheme and by choosing the renormalization 
scale \i = A = m K , we obtain the following coefficients 



< = ^f^, (25) 



< = ^|^. (26) 

Adding up the contributions from these operators and the contributions com- 
ing from the diagram in fig. 1, i.e., substituting eqs. (25), (26) in eqs. (17)-(20), 
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K ++ 

c 

e kR 



-®- 



e iR 



e iR 



K ++ 

^--^ 

c 

e kR 



(a) 



"jR 



e iR 



c 

e kR 



iR »» 



K 



e iR 



-jR 



e iR 



(b) 



Fig. 3. One-loop matching. 

we obtain exactly the same amplitudes as obtained from a calculation in the 
full model up to terms C((m M /m K ) 4 ). Therefore the full amplitude for /i-e 
conversion is dominated by the divergent contribution from the diagram in 
fig. 1. In the effective field theory language one says that the amplitude is 
dominated by the running from the scale of new physics A = m K to relevant 
scale of the process: in the case of /i-e conversion, m M (or m r if r leptons are 
running in the loop). This conclusion is independent on the model as long as 
four-fermion interactions eq. (6) and/or eq. (7) exist. 



As expected, in this model there are no contributions to e^-y^e^-Z^ since the 
physics of a scalar singlet decouples for m K ^> A F . This means that the con- 
tributions of the scalar can only lead to operators suppressed by, at least, 
For instance one could easily obtain an operator like eq. (3) with the 
photon replaced by the Z-boson. The contribution of those operators to /i-e 
conversion, however, are suppressed by a factor (m^/m z ) 2 with respect to the 
photonic contributions. 
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5 Numerical results and conclusions 

Now we are ready to compare the branching ratios of /i-e conversion and 
fi — > e 7 in the class of models we consider and to analyze the relative potential 
of different experiments to test muon flavour conservation. For definiteness we 
consider only the right-handed operators in lagrangian (1). We first study 
the general case in the framework of our effective theory and then we apply 
the results to our specific model. Substituting the form factors (17)- (20) to 
eq. (14) and taking only the dominant logarithmic terms (with A = 1 TeV in 
the logarithms) we find 



where the first number in the expressions corresponds to k — e, fi and the 
number in the brackets to k — r. If k — r the logarithm enhancement is 
m(m T /A) instead of ln(m M /A) and it is slightly smaller than in the \x or e 
cases. Note that the conversion is somewhat enhanced in Pb and Au (in fact, 
maximized [25]) if compared with T%. 

In the effective lagrangian framework /i — > e 7 does not get contributions from 
loops in fig. 1 and, therefore, it is not enhanced by large logarithms. In any 
full theory in which both /i-e conversion and ji — > e 7 are induced by loops all 
the couplings should be of the same magnitude (compare, e.g., eq. (24) with 
eq. (26) in our doubly charged scalar model). Assuming a aR = a RR we obtain 



for any type of fermion in the loop. 

Comparison of eqs. (27), (28) with eq. (29) shows that due to the presence 
of large logarithms the [i-e conversion rate is comparable or even exceeds the 
H — > e 7 rate. To constrain new physics we have to also take into account 
the sensitivity of experiments. The present experimental upper limits on the 
branching ratios of the processes are Rj^(exp) ^ 4.3 • 1CT 12 [1], R P ^(exp) ^ 
4.6 • KT 11 [2] and R^exp) £ 4.9 • KT 11 [31]. SINDRUM II experiment at PSI 
taking presently data on gold will reach the sensitivity pj^,ex-pected < 5 . io~ 13 
[30] and starting next year the final run on Ti it should reach R^i, expected < 
3- 10~ 14 [4]. Normalizing the branching ratios to the experimental upper limits 
we get for A = 1 TeV 




(27) 



(28) 




(29) 
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R = 


4.6- 10~ n 


4.3- i(r 12 


5.0- 1(T 13 


3.0- 10~ 14 


log-enhanced a-e 


32 


44 


101 


158 


non-enhanced a-e 


7 


9 


20 


32 




23 


41 


70 


141 



Table 1 



Values of A, in TeV, probed in u-e conversion and a — ► e 7 for different upper bounds 
on the branching ratios. The upperbound 4.6 • 10 -11 is the present bound for a-e 
conversion on Pb and it is very close to the present /j^e] bound (4.9 • 10~ n ). 
4.3 • 10~ 12 is the present bound for u—e conversion on Ti. 5 • 10~ 13 and 3 • 10 -14 are 
the expected bounds in the next year for u-e conversion on Au and Ti, respectively. 



R 



Ti 
fie 



4.3- 10- 12 

p^Pb, Au 



"fie 



4.6-10 



-11 



11.4 (5.7) B Ti 
3.1 (1.25) B Au 



R~, 



4.9 • 10- 11 



4.9-10 



-11 



(30) 
(31) 



where, again, numbers in the brackets correspond to the case k = r and the 
factors B, B = R 1 ^ e esent (exp)/Rj l ^ ture (exp), take into account the improve- 
ments in the experimental sensitivity. Eq. (30) and eq. (31) constitute the 
central result of this work: in the class of models we consider searches for 11-e 
conversion in both Ti and Pb constrain new physics more stringently than 
searches for \x — > e 7. In addition, from SINDRUM II one expects to achieve 
B Au = 92. already in forthcoming months and B T% = 1.4- 10 2 next year. If the 
aimed sensitivity will be achieved then \i-e experiments probe the couplings 
a of new physics more than one order of magnitude more stringently than 
fi -> ery. 



To show which scales of new physics A can be probed in \i—e conversion and 
/i^e) experiments we have presented the values of A in TeV-s in Table 1 for 
different experimental upper bounds on the branching ratios of the processes. 
All the couplings a are taken to be equal to unity. We have considered both 
classes of models with and without logarithmic enhancement of fx-e conversion. 
If the experimental limits for 11-e conversion and 11 — > e 7 are equal then 11-e 
conversion enhanced by large logarithms has better sensitivity to A than 11 — > 
e 7, especially in the case of Pb and Au experiments. The scales testable reach 
A ~ O(10 2 ) TeV. However, fi-e conversion without logarithmic enhancement 
can only probe scales lower by about a factor 5. 



To illustrate the discussion above let us present the experimental bounds on 
the couplings h of k ++ in our model. Substituting the couplings in eqs. (24)- 
(26) to eqs. (17)- (20) and using the present experimental limit for Ti we obtain 
from Li-e conversion for m K = 1 TeV 



12 



6 • 1(T 4 



e/i'^ee ' MM e/z ~ 



9 • 1(T 4 




(32) 



while /i — > e7 gives 



^k^e* £3-10 



,-3 



(33) 



The bounds (32) are new limits on the off-diagonal doubly charged scalar 
interactions (note that tree level \x — > 3e probes only h^h*^). While derived 
for the right-handed singlet the limits apply with a good accuracy also for the 
interactions of triplet scalars appearing in models with enlarged Higgs sectors 
as well as in left-right symmetric models This is because the doubly charged 
component of triplet gives the dominant contribution both to //-e conversion 
and /i — > e 7. Note that the upper bounds (32) are going to be improved by 
an order of magnitude with new /i-e conversion data. 

Finally, we would like to stress that our main result, the logarithmic enhance- 
ment of /i-e conversion rate, is completely general and applies to all models 
with effective interactions of four charged fermions. For simplicity we have 
constrained ourselves to purely leptonic operators. However, the same effect 
is also present for operators involving quarks. To get large logarithms one just 
needs light charged fermions in the loop. Therefore, loop induced /i-e conver- 
sion is also enhanced in models with broken _R-parity [32] and leptoquarks but 
not in _R-conserving MSSM or SUSY GUT's considered in ref. [7] in which the 
light fermions in loops are necessarily neutral. 

In conclusion, using the effective lagrangian description of new physics we have 
pointed out a wide class of models with effective four charged fermion inter- 
actions in which loop induced /i-e conversion in nuclei is enhanced by large 
logarithms. With the present upper limits on /i-e conversion and /j, — > 
branching ratios bounds on new physics (occurring at loop level) derived from 
these processes are more restrictive in the case of /i-e conversion. In nearest 
future this factor will increase by more than one order of magnitude due to the 
expected improvements in sensitivity of already running /i-e conversion exper- 
iments. This general result is confirmed by exact calculations in the extension 
of the SM with doubly charged singlet scalar. 
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